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ABSTRACT 


Theoretical performance, both equilibrium and rate-limited, is presented for the 
oxygen difluoride -diborane propellant system. A one -dimensional kinetic flow model 
was used. The work also includes identification of the rate controlling chemical reac- 
tions in the recombination process. Reactions are quantitatively rated in importance 
and the sensitivity of performance to the uncertainty of individual rate constants is 
studied in detail. Results indicate that the hydrogen-atom recombination is the pre- 
dominant reaction in the entire process as far as performance is concerned and that no 
reaction involving boron species has a significant effect on performance. The theoret- 
ical performance level of oxygen difluoride-diborane is significantly higher than that for 
the methane- fluorine /oxygen propellant, the most energetic of the hydro carbon- 
fluorine /oxygen systems. 



THEORETICAL NONEQUILIBRIUM PERFORMANCE OF OXYGEN 
DIFLUORIDE-DIBORANE ROCKET PROPELLANT 
by David A. Bittker 
Lewis Research Center 

SUMMARY 

Theoretical computations have been performed for the oxygen difluoride-diborane 
propellant system to obtain both equilibrium and rate-limited theoretical performance 
over an oxidant-fuel ratio range of 2. 7 to 4 at chamber pressures of 100 and 300 pounds 
per square inch absolute (68. 95 and 206. 9 N/sq cm). The performance limiting reac- 
tions in the recombination process were identified and quantitatively rated in importance 
A detailed study was also made to determine the sensitivity of vacuum specific impulse 
to individual uncertainties in the individual rate constant data. Finally, the theoretical 
performance of oxygen difluoride was compared with that for two hydrocarbon-fluorine/ 
oxygen systems. 

Results of these computations showed that kinetically limited vacuum specific im- 
pulse at 100 psia (68. 95 N/sq cm) chamber pressure ranged from 411 to 422 pound- 
seconds per pound (4030 to 4138 (N-sec)/kg) for oxygen difluoride-diborane. This is 
considerably higher than the maximum value 398 pound-seconds per pound (3903 
(N-sec)/kg) for the methane-fluorine/oxygen system at the same chamber pressure. The 
most important reaction occurring in the nozzle is the hydrogen atom recombination 
process. Performance is very sensitive to uncertainty in the rate constant for this re- 
action over the entire O/F range. To a much smaller extent performance is also af- 
fected by uncertainties in the rate constant for the H + F and H + OH recombinations. No 
reaction involving boron species is rate controlling. In fact all such reactions can have 
at least one order of magnitude uncertainties in their rate constants without causing any 
change in the computed performance values. These results are in agreement with the 
predictions of the reaction rating parameter X H •, which gives both the magnitude and 
direction of the enthalpy change rate caused by a particular reaction. 


INTRODUCTION 


In the continuing development of new chemical rocket propellant systems interest 
has centered on fluorine containing combinations other than the hydrogen-fluorine system 
itself. For many missions a less energetic propellant than Hg-Fg, but one with fewer 
handling and storage problems, can be used. Several of these systems are termed 
"space storables, " inasmuch as the constituents can be stored indefinitely as liquids in 
outer space with use of only insulating materials and heat shields. Examples of these 
space storable fuels are the light hydrocarbons (e.g. , methane, propane, and butene- 1) 
and also diborane (BgHg). These fuels are usually combined with fluorine /oxygen mix- 
tures or oxygen difluoride (OFg) as an oxidizer. A detailed study of the nonequilibrium 
performance of two hydrocarbon-fluorine/oxygen systems is given in reference 1. This 
work also presents equilibrium performance curves which compare several space- 
storable propellants with the best earth storable and cryogenic systems in use at the 


System Chamber pressure, 
psia (N/cn/l 



Figure 1. - Performance of rxket propellants assuming equilibrium expansion to area ratio 
of 40. 
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present time. For reference purposes some of these curves are reproduced in figure 1. 
A similar comparison is shown in a recent NASA survey paper (ref. 16) which also gives 
a comprehensive discussion of the NASA space storable propulsion programs, including 
objectives and future goals. The hydrogen-fluorine and Aerozine- 50-nitrogen tetroxide 
combinations are used as the current upper and lower performance limits. The space 
storable fuels shown are diborane, combined with oxygen difluoride and with oxygen/ 
fluorine mixtures, and the hydrocarbons methane and butene- 1 combined with oxygen/ 
fluorine mixtures. These last two systems are the ones studied in reference 1. These 
equilibrium results show that diborane, with either OFg or oxygen/fluorine mixtures, 
has a higher performance level than methane, the most energetic hydrocarbon fuel. 
However, the best space storable system still has considerably less impulse than 
hydrogen-fluorine. 

A study of the OFg-diborane system was undertaken because of its favorable per- 
formance and because of renewed interest in it as a practical rocket propellant. The 
main effort was the computation of rate-limited performance and the identification of the 
important chemical reactions occurring during the nozzle expansion process. A method 
of quantitatively rating reactions in order of importance was given in reference 1 and has 
been used here. In addition we have studied the effect on computed specific impulse of 
changing individual rate constants over the range of their estimated uncertainty. A com- 
plete set of equilibrium, kinetic and frozen performance computations is presented for a 
wide range of operating conditions for the OFg-diborane system. The calculations were 
made for chamber pressures of 100 and 300 psi (68. 95 and 206. 9 N/ sq cm) and an 
oxidant-fuel ratio (O/F) range of 2. 7 to 4. A comparison is also given of the rate- 
limited performance for OFg-diborane with that for methane-fluorine/oxygen and 
butene- 1 fluorine oxygen. 


PROCEDURE 

Equilibrium and Frozen Computations 

Equilibrium and frozen computations were performed to obtain the starting condi- 
tions for the kinetically- limited computations. The necessary conditions are the signif- 
icant species present at equilibrium in the convergent nozzle section, their concentra- 
tions, and the equilibrium values of the fluid dynamic variables, velocity, temperature, 
and either pressure or density. From these computations we also obtain the maximum 
and minimum values of theoretical performance for comparison with the rate-limited 
performance. The significant species were found by omitting several species which pre- 
liminary equilibrium calculations showed to be present with a mole fraction of less than 
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0. 001. These trace species were dropped from consideration after it was verified that 
these omissions caused negligible change in all equilibrium performance results up to 
an expansion area ratio of 40 in the nozzle. A negligible change in vacuum specific im- 
pulse was arbitrarily taken as ±0. 5 pound-second per pound (±4. 9 (N-sec)/kg) or less. 
All equilibrium and frozen performance results presented in this report were computed 
using the method and program of Zeleznik and Gordon (ref. 3). 

Kinetical ly-Limited Computations 

After selecting the significant species present in the gas mixture entering the 
nozzle, a reaction mechanism was hypothesized for the complete nozzle recombination 
process. Then a set of "reference” rate constant expressions was selected by survey- 
ing the available literature. Many of these reactions are well known ones for which ex- 
perimental rate constant data are available. Practically no experimental work has, 
however, been done on reactions of the boron containing species. Therefore, it was 
necessary not only to hypothesize the reactions which could occur, but also to calculate 
approximate theoretical rate constant expressions for these reactions. The details of 
the entire reaction and rate constant selection procedure are given in a subsequent sec- 
tion. Computations were made using these reference rates for all assigned experi- 
mental conditions to obtain a set of reference performance curves. These results also 
showed which reactions are important in determining kinetic performance. A detailed 
rate variation study was then done to determine the sensitivity of performance to indi- 
vidual rate constant variations. This type of study is especially important for the OFg- 
diborane system since much of the rate constant data have large uncertainties. 


THEORY AND INPUT DATA 
Kinetic Flow Model and Equations 

The problem of calculating chemical-rate limited rocket performance is just a par- 
ticular example of the general problem of following the progress of several simultaneous 
chemical reactions in a flow system. One must solve a series of simultaneous differen- 
tial equations consisting of fluid dynamic conservation equations coupled with the equa- 
tions for the rates of change of species concentrations due to chemical reaction. The 
particular set of equations used for this work is based on the one dimensional flow of a 
uniform mixture of thermally perfect gases resulting from complete combustion of 
rocket propellants to chemical equilibrium in the chamber. All losses due to viscosity 
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and heat transfer are neglected. The specific system of equations used was first given 
in reference 2. It consists of the following N + 3 equations in the N + 3 unknowns V, T, 
and p, and the concentrations of the N chemical species. All symbols are defined in 
the appendix. 


dY = ..V. (jdA-j, 

dx jfl _ ^ \A dx 


( 1 ) 


dT 

dx 


= -T 


(r - DJr 




2 . j \A dx 


( 2 ) 


dx 


Jt 2 

Jl 2 - 1 


\A dx j 


+ si 


(3) 


da. w.r* 

— = ^— i = 1, 2, . . ., N (4) 

dx P V 


The auxiliary equations needed are 


p = ^T 
M 


(5) 


Jt 


2 _ V*M 
yRT 


( 6 ) 


r = 


R 

M 
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The si and terms are given by 


N 


si = 


r*T 

PV 


YJs'- 


( 8 ) 


i=l 
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Vi 


Equations (1) to (4) constitute a system of differential equations of the form 


— = f i( x , Yj, y 2 . ■ ■ o y n ) i = 1, 2, . . n (10) 

The method of solution used for this system is the implicit numerical integration 
technique described by Tyson and Kliegel (ref. 4). This technique was used in a per- 
formance program written for NASA Manned Spacecraft Center. The present computa- 
tions were performed with an extensively modified version of that program. As is cus- 
tomary in all flow-kinetic programs, chemical reactions are considered reversible and 
the law of microscopic reversibility is used to calculate the reverse reaction rate con- 
stant from the forward rate constant and the equilibrium constant for the reaction. 


Reaction Mechanism and Species 


Equilibrium flow computations showed that the following significant species are 
present in the expanding gas mixture for the OFg-diborane propellant system. 

B, O, H, F, BO, H 2 , 0 2 , HF, H 2 0, OH, BOF, BF, B0 2 , BF 2 


and BFg. 

The reactions selected for the complete reaction mechanism for the recombination 
process are listed in table I along with their heats of reaction at 298 K. 

The reference rate constants selected for these reactions are given in table n. 
Those for reactions I to IX are identical to the rate constants used for the hydrocarbon- 
F 2 /0 2 study reported in reference 1. These rate constant expressions are selected val- 
ues based on an evaluation of all the available data. The rate constant expressions for 
the second group, which includes all other reactions, are theoretical estimates usually 
taken either from reference 5 or 14. The rate constants for reactions XI and XII, how- 
ever, were calculated during the present work using the modified collision theory of 
C. Fredrick Hansen as detailed in reference 6. The accuracy of this second group of 
rate constants must be considered order of magnitude only. Therefore, the study of the 
effect of rate constant variation on performance is especially important for these reac- 
tions. 
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TABLE I. - CHEMICAL REACTIONS IN THE NOZZLE RECOMBINATION 


PROCESS FOR THE OXYGEN DIFLUORIDE- DIBORANE 


PROPELLANT SYSTEM 


Reaction 

number 

Reaction 

Heat of reaction 
(forward direction), 

AH 298 

kcal/mole 

kJ/mole 

I 

H+F + M-HF + M 

-136 

-569 

n 

h + h + m-h 2 + m 

-104 

-435 

in 

F + H 2 - HF + H 

-32 

-134 

IV 

H + OH + M “ HgO + M 

-119 

-498 

V 

O+H + M-OH + M 

-102 

-427 

VI 

H 2 + OH — H 2 0 + H 

-15 

-63 

VII 

H + 0 2 ~ OH + O 

+ 17 

+71 

vra 

O + H 2 ~ OH + H 

+2 

+ 8 

rx 

o + o + m-o 2 + m 

-119 

-498 

X 

F + OH Z. O + HF 

-33 

-138 

XI 

BO + F + M ~ BOF + M 

-169 

-707 

XII 

BF + O + M ~ BOF + M 

-160 

-669 

XIII 

BO + HF ~ BOF + H 

-33 

-138 

XIV 

BF + OH U BOF + H 

-58 

-243 

XV 

B + O + M-BO + M 

-189 

-791 

XVI 

B+F + M-BF + M 

-197 

-824 

XVII 

BF 2 + BO - BOF + BF 

-58 

-243 

xvin 

BF 2 + H ~ HF + BF 

-25 

-105 

xtx 

O + BF ~ BO + F 

+ 8 

+33 

XX 

BF + B0 2 ~ BOF + BO 

-10 

-42 

XXI 

BF 2 + HF ~ BF 3 + H 

-18 

-75 

XXII 

BF + F + M — BF 2 + M 

-111 

-464 

xxm 

BF 2 + F + M = BFg + M 

-154 

-644 

XXIV 

BO + O + M - B0 2 + M 

-150 

-628 

XXV 

BF + Og — BOF + O 

-41 

-172 


TABLE n. - REFERENCE REACTION RATE CONSTANTS 3 -’ b 


Reaction 

Reaction 

Preexponential 

Temperature 

Activation 

Reference 

number 


factor, 

exponent, 

energy, 




D 

n 

E a’ 

cal /mole 


I 

H+F + M-HF + M 

3. 0xI0 1S 

-1. 0 

0 

7 

n 

h + h + m~h 2 + m 

7. 5xl0 18 

-1.0 

0 

7 

m 

F + H 2 ~ HF + H 

2.04X10 13 

-. 12 

3 750 

7 

IV 

H + OH + M ~ H 2 0 + M 

4. 5X10 21 

-1. 5 

0 

8 

V 

0 + H + M~0H + M 

4. OxlO 18 

-1.0 

0 

9 

VI 

vn 

H 2 + OH ~ H 2 0 + H 
H + 0 2 = OH + O 

2. 3xl0 13 
2. 2XI0 14 

0 

0 

5 200 
16 600 

10 

11 

vra 

O + H 2 - OH + H 

4. OxlO 13 

0 

10 200 

12 

IX 

o + o + m-o 2 + m 

8. 15X10 18 

-1. 22 

0 

13 

X 

F + OH - O + HF 

2. 9X10 12 

. 68 

200 

14 

XI 

BO + F + M ~ BOF + M 

2.5X10 18 

-.77 

-843 

(c) 

xn 

BF + O + M : BOF + M 

4.77X10 18 

-. 87 

-443 

(c) 

xin 

BO + HF - BOF + H 

1. OxlO 11 

. 5 

7 000 

5 

XIV 

BF + OH “ BOF + H 

1. OxlO 11 

. 5 

6 000 

5 

XV 

B+O + M-BO + M 

3. OxlO 16 

-.5 

0 

5 

XVI 

B + F + M — BF + M 

3. OxlO 16 

-. 5 

0 

5 

XVII 

BF 2 + BO “ BOF + BF 

1. OxlO 11 

.5 

6 000 

5 

xvra 

BF 2 + H — HF + BF 

1. 05X10 11 

. 5 

6 663 

14 

XIX 

O + BF - BO + F 

5. 3X10 11 

. 5 

9 865 

14 

XX 

BF + BOg - BOF + BO 

1. 05X10 11 

. 5 

7 585 

14 

XXI 

XXII 

BFg + HF “ BFg + H 
BF + F + M — BF„ + M 

1.05X10 11 
1. OxlO 11 

. 5 
. 5 

7 468 
0 

14 

5 

xxm 

XXIV 

i 

BFg + F + M ~ BFg + M 
BO + O + M - BOg + M 

I. OxlO 15 
1. OxlO 15 

-. 5 
-. 5 

0 

0 

14 

14 

XXV 

BF + Og - BOF + O 

6. 96X10 11 

. 5 

0 

14 


a Rate constant expression is k = DT n exp(-E„/RT). 

D Units of k are cm J /(mole)(sec) for a bimolecular reaction and cm D /(mole)‘ s (sec) 


for a three- body reaction. 
c Computed using the theory of ref. 6. 



RESULTS AND DISCUSSION 


Performance Results 

Theoretical performance calculations were made for chamber pressures of 100 and 
300 psi (68. 95 and 206. 9 N/sq cm) over the oxidant-fuel ratio range of 2. 7 to 4, the ra- 
tios of greatest interest to engine designers. The stoichiometric equation for the com- 
plete consumption of all reactants is 

30F 2 + B 2 H 0 - 4HF + 2BOF + HgO 

The stoichiometric O/F is 5. 851 and therefore, engines using this propellant system 
operate in the fuel-rich region. Computed performance for the equilibrium, kinetically 
limited and frozen assumptions is given in figure 2 for the two chamber pressures. A 
typical 15° conical nozzle with throat area of 8. 65 square inches (55. 81 sq cm) and an 
expansion ratio of 40 was used for the kinetically limited computations. Thrust level is 
approximately 1500 pounds at 100 psi chamber pressure (6672 N at 58.95 N/sq cm). All 
computations are based on the latest experimental value of the heat of formation of oxy- 
gen difluoride (5. 86 kcal/mole or 24. 52 kJ/mole) reported by King and Armstrong in ref- 
erence 15. Theoretical performance is moderately affected by changing O/F. Figure 2 



Oxidant/fuel weight ratio, 0/F 

Figure 2. - Theoretical performance of oxygen difluoride-diborane. 
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shows that the kinetically limited performance increases steadily as O/F increases. 
Specific impulse ranges from 411 to 422 pound- seconds per pound (4030 to 4138 
(N-sec)/kg) at a chamber pressure of 100 psia (68. 95 N/sq cm). At the higher chamber 
pressure, the nonequilibrium specific impulse 1^ is between 423 and 434 pound-seconds 
per pound (4148 and 4256 (N-sec)/kg). Increasing chamber pressure, therefore, in- 
creases 1^ about 12 units or 3 percent. The kinetic impulse loss I - 1^ is reduced 
from about 16 units (at 100 psia) to approximately 8 units by raising the chamber pres- 
sure. The increased molecular collision rate significantly speeds up the recombination 
reaction rates at the higher pressure. 

The dependence of kinetic loss on position in the nozzle is shown in figure 3, where 
the three types of theoretical performance are plotted as a function of supersonic area 
ratio for a typical expansion. The assigned conditions are a chamber pressure of 100 psi 



Figure 3. - Theoretical performance of oxygen difluoride-diborane plotted against area ratio. 
Chamber pressure, 100 pounds per square inch absolute (68. 95 N/cm 2 ); oxidant-fuel 
ratio, 3. 6. 
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and an oxidant-fuel ratio of 3. 6. The figure shows that kinetic specific impulse deviates 
from the equilibrium value by no more than one impulse unit up to an area ratio of 3. 

The deviation is 4 units or 1 percent of the equilibrium value at an area ratio of about 6. 
The kinetic loss at the final area ratio of 40 is 14 units or about 3 percent of the equilib- 
rium specific impulse. In reference 1 a typical methane-fluorine/oxygen expansion at 
the same chamber pressure showed a 1 percent kinetic loss at an area ratio of 3. 5 and 
a total loss of 5 percent at area ratio 40. This single comparison, using the same 
chamber pressure, but different nozzles, indicates that the OFg-diborane combination 
has significantly less kinetic loss than the methane-fluorine/oxygen system. 


Effect of Reaction Rate Variation on Performance 


A general idea about the relative importance of the many chemical reactions used in 
the assumed mechanism can be obtained from the changes in gas composition. Table III 
shows species mole fractions in the nozzle gases at the throat and at the end of the 

TABLE ra. - NOZZLE GAS COMPOSITION FOR THE 
OXYGEN DIFLUORIDE-DIBORANE PROPELLANT 


[Chamber pressure, 100 psia (68. 95 N/sq cm). ] 


Species 


Oxidant fuel ratio 




2.7 

40 



Species mole fraction 


Throat 

Area ratio, 40 

Throat 

Area ratio, 40 

B 

0.00101 

0.00106 

0. 00036 

0.00038 

BF 

. 09663 

. 10183 

.02531 

.00398 

BF 2 

.00656 

.00032 

.00355 

3. 6X10" 5 

bf 3 

.00168 

.00368 

.00105 

. 00031 

BO 

.01940 

.00019 

.02115 

. 00031 

B °2 

.00126 

1. 5xl0 -6 

. 00898 

.00369 

BOF 

. 20472 

.23944 

. 21025 

. 27862 

F 

.00134 

4. lxlO" 7 

.01951 

6. 4X10" 5 

H 

. 20224 

0. 14063 

. 28222 

. 25478 

HF 

. 13774 

. 12673 

.28968 

. 30551 

H 2 

. 32477 

. 38589 

. 10986 

. 14382 

h 2 ° 

o 

.00189 

.00021 

.00023 

4.4X10" 7 

.00576 

.01306 

. 00680 
.00106 

OH 

. 00057 

1. lxlO" 6 

. 00898 

. 00061 

°2 

4.7X10' 7 

1.5X10' 10 

. 00027 

1.9X10" 5 
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nozzle. For O/F = 2.7, the major species are BF, BOF, H, HF, and H^. The signifi- 
cant changes among these species are the increase in BOF and Hg mole fractions and the 
decrease in H and HF mole fractions. At O/F = 4, the major species are BOF, H, HF, 
and Hg. The same changes in H, BOF, and Hg concentrations occur at O/F = 4 as did at 
O/F = 2. 7. But HF concentration now increases between throat and nozzle exit instead 
of decreasing. So one would look for reactions that form E^ and BOF and consume H as 
being the most rate controlling processes; however, a more quantitative and systematic 
procedure is needed to identify these reactions. The method to be used here was de- 
scribed in reference 1 and is a modification of the technique originally proposed in ref- 
erence 14. A reaction rating parameter Xjj . is used to quantitatively measure the 
importance of each individual reaction. It is proportional to the extent the reaction pro- 
ceeds (i. e. , the net volumetric conversion rate) and to the heat of reaction. The ex- 
pression for Xtt i is 


X H, j - X j^ AH 298^ 


( 11 ) 


where 


X j 


is a net reaction conversion rate parameter given by the equation 


X . 
X. =-£J 
1 2 
P 



P 


2 


( 12 ) 


The quantity X . is the net volumetric conversion rate for the general reaction: 
^5 J 


N 


y u..s. 

Lj lj i 
i=l k 


k fj N 


bj 


y u..s. 

Lj i] i 
i=l 


(13) 


4 th 


where S. represents the name of the i species in the system. The rating parameter 

Xjj . measures the net enthalpy change rate caused by a particular reaction. A negative 

Xtt , means that energy is released into the gas stream, while a positive value means 
J 

that energy is absorbed due to the reaction. Calculated values of X n . for the reac- 

H, J 

tions under consideration are shown in table IV for two O/F values and a chamber pres- 
sure of 100 psia (68. 95 N/sq cm). At O/F = 2.7 the hydrogen atom recombination ap- 
pears to be much more important than any other reaction, since it releases the largest 
amount of energy into the reacting gases. This recombination reaction is also the most 
important one at the high O/F value of 4.0. But several other reactions now appear to 
be quite important in releasing energy. It should be noted that the H + OH recombination 
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TABLE IV. - QUANTITATIVE RANKING OF REACTIONS FOR OF„-B„H„ NOZZLE 

1 u 


EXPANSION USING THE ENTHALPY EXCHANGE RATE PARAMETER X„ . 
[Chamber pressure, 100 psia (68. 95 N/sq cm). ] 


Rank 

Oxidant fuel ratio 


2.7 

4 


Reaction 

X H,i 

r h,J 

Reaction 

X H,J 

R H,j 

1 

h + h + m-h 2 + m 

-403 800 

1. 0000 

h+h + m-h 2 + m 

-311 300 

1.0000 

2 

BO + HF - BOF + H 

-60 850 

. 1507 

H+F + M- HF + M 

-117 800 

. 3784 

3 

H+F + M-HF + M 

-14 130 

.0350 

H + OH + M - H 2 0 + M 

-101 900 

. 3273 

4 

H + OH + M - H 2 0 + M 

-12 410 

.0307 

BO + HF - BOF + H 

-100 800 

. 3238 

5 

BF + OH — BOF + H 

-10 120 

.0251 

BF + OH - BOF + H 

-24 670 

.0792 

6 

H 2 + OH - HgO + H 

+6 666 

.0165 

F + H 2 - HF + H 

-23 320 

.0749 

7 

BO + BF 2 - BOF + BF 

-6 328 

.0157 

H 2 + OH - H z O + H 

+ 16 200 

.0520 

8 

F + OH ~ HF + O 

-2 229 

.0055 

BF + O — BO + F 

+ 11 370 

.0365 

9 

H + BF 2 - HF + BF 

-1 971 

.0050 

H + BF 2 ~ HF + BF 

-8 423 

.0271 

10 

B0 2 + BF “ BOF + BO 

-1 835 

.0045 

O+H + M-OH + M 

-7 889 

.0253 

11 

BO + F + M - BOF + M 

-1 131 

.0028 

BO + F + M - BOF + M 

-7 576 

.0243 

12 

BF 2 + HF — BFg + H 

-1 089 

.0027 

BF + 0 2 ~ BOF + O 

-6 183 

.0199 

13 

BF + 07 BO + F 

+ 1 022 

.0025 

BF + O + M — BOF + M 

-5 528 

.0178 

14 

F + H 2 - HF + H 

-1 002 

.0025 

F + OH — HF + O 

+5 283 

.0170 

15 

BF + O + M - BOF + M 

-789 

.0020 

BO + BF 2 - BOF + BF 

-3 994 

.0128 

16 

O + H + M — OH + M 

-212 

.0005 

B0 2 + BF - BOF + BO 

-3 399 

.0109 

17 

O + H 2 - OH + H 

-132 

.0003 

O + H 2 " OH + H 

+ 2 350 

.0075 

18 

BF + 0 2 - BOF + O 

+48 

.0001 

h + o 2 r OH + O 

-2 120 

.0068 

19 

H + 0 2 - OH + O 

-32 

.0001 

BF 2 + HF - BFg + H 

+ 1 946 

.0063 


Relative rate for H atom 



Oxidant/fuel ratio, 0/F 


Figure.4. - Effect of rate ofH+H + M-H 2 + M reaction on 
performance of oxygen dif I uor ide-dibora ne. Chamber 
pressure, 100 pounds per square inch absolute 
(68.95 N/cm^f; final expansion ratio, 40. 



reaction has a very high Xjj . value even though OH and HgO are minor species. Its 
high specific rate constant and heat of reaction counterbalance the effect of the small 
concentrations to make its net energy release effect quite large. 

The exact effect of each important reaction was determined by individually varying 
its specific rate constant over the range of its estimated uncertainty. The results of 
this variation for the H atom recombination are shown in figure 4. The value of kg was 
changed from five times to one fourth of the reference value. It is clear that variation 
of this rate constant has a large effect on specific impulse over the entire range of as- 
signed conditions used. Changing the rate constant by its uncertainty changes perform- 
ance by 3 to 6 pound-seconds per pound (29 to 58 (N-sec)/kg). The effect of changing 
other reaction rate constants is summarized in table V. Computations were made only 
at the highest and lowest O/F values for a chamber pressure of 100 psia (68. 95 N/sq cm). 
The main conclusion from table V is that only two rate constants other than the one for 
the H + H reaction have a large enough uncertainty to effect performance by more than 
0. 5 pound-seconds per pound (4.9 (N-see)/kg). These are the constants for 
H + F + M ~ HF + M and H + OH + M ~ HgO + M. This effect is only apparent in the 
high O/F region. Another fact brought out in these rate variation studies is that no re- 
action among boron- containing species is really rate controlling on the entire process. 
The two most important ones (BO + HF ~ BOF + H and BF + OH — BOF + H) are low en- 
ergy shuffle reactions, whose rate constant variation by one order of magnitude causes 
negligible change in performance. The two most energetic reactions involving the BOF 
species are not even among the first 10 significant reactions. These are the three-body 
recombinations of BO and F or BF and O to give BOF. A few rate variation computa- 
tions were, however, made for these reactions to determine the effect if their rate con- 
stants were in error by two orders of magnitude. The results, shown in table V, indi- 
cate that the effect of both reactions would still be insignificant under these conditions. 

The last column of table V shows the effect on performance of removing each reac- 
tion individually from the mechanism. Only for one reaction (BO + HF — BOF + H) does 
this removal cause a significant change in performance from the reference value. This 
fact does not mean, of course, that several of these reactions can be removed at the 
same time without changing our results. As in reference 1, we shall compare the reac- 
tions in a given mechanism with one another by using a ratio of Xtt i values, namely: 




(14) 


In this equation (Xtt j is the largest value for all reactions considered. The 
v ft, J/ max 

smaller the R H • value, the less important a reaction should be. Reference 1 shows 
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TABLE V. - EFFECT OF REACTION RATE VARIATION ON PERFORMANCE OF 


OXYGEN DIFLUORIDE-DIBORANE 


[Chamber pressure, 100 psia (68.95 N/sq cm). ] 


Number 

Reaction 

Uncertainty factor 

Oxidant 

fuel 

ratio, 

O/F 


Rate constant 


Reference 

Vai 

Maximum 3, 

Minimum* 3 

0 


Maximum 

Minimum 

cuum specific impulse 

I 

H + F + M — HF + M 

4 

0.25 

2.7 

411.4 

411. 5 

411.4 

411.4 





4 

422.0 

422.7 

421.9 

421. 8 

IV 

H + OH + M - H 0 0 + M 

10 

.10 

2.7 

411.4 

411.7 

411.4 

411.4 





4 

422.0 

424.7 

421.6 

421.5 

xra 

BO + HF - BOF + H 

10 

. 10 

2.7 

411.4 

411. 5 

411. 1 

410.4 





'4 

422.0 

422. 1 

421.7 

420. 1 

XIV 

BF + OH ~ BOF + H 

10 

. 10 

2.7 

411.4 

411. 5 

411.4 

411.4 





4 

422.0 

422. 3 

422.0 

422.0 

VI 

H„ + OH - H„0 + H 

2 

.5 

2.7 

411.4 

411.4 

411.4 

411.2 





4 

422.0 

422. 1 

422. 1 

421. 8 

m 

F + Hg — HF + H 

10 

. 10 

2.7 

411.4 

411.4 

411.4 

411.4 





4 

422.0 

422.0 

422. 1 

422. 1 

X 

F + OH - HF + O 

10 

. 10 

2.7 

411.4 

411. 5 

411.5 

411. 5 





4 

422.0 

421.7 

422. 1 

422. 1 

XIX 

BF + O - BO + F 

10 

. 10 

2.7 

411.4 

411.5 

411.4 

411.4 





4 

422.0 

422.4 

421.9 

421. 8 

XI 

BO + F + M - BOF + M 

100 

.01 

2.7 

411.4 

411.6 

411.4 






4 

422.0 

422. 5 

422. 1 


XII 

BF + O + M - BOF + M 

100 

.01 

2.7 

411.4 

411. 5 

411.4 






4 

422.0 

422. 5 

422. 1 



a Reference rate constant x maximum uncertainty factor. 
^Reference rate constant x minimum uncertainty factor. 


that all reactions with an FUr . of about 0. 02 or less can be removed as a group from 

n , 3 

the hydrocarbon-fluorine /oxygen mechanism without changing kinetic specific impulse 

by more than the arbitrary limit of ±0. 5 pound-second per pound. Reference 1 also 

shows that an individual reaction with an R™ ■ of 0. 14 can be removed without changing 

1 

performance more than this arbitrary limit. For the OFg-diborane system the Rjj j 
values in table IV and the information in table V show the following results. At 
O/F = 2.7, any reaction with R^ ^ less than or equal to 0.035 can be individually re- 
moved without significant change in performance. For O/F = 4, the reaction 
H + F + M ~ HF + M with R H , equal to 0. 3782 can be removed with negligible change. 
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However, when the reaction BO + HF — BOF + H with Rtt ■ equal to 0. 3238 is removed, 

J 

specific impulse is decreased by 1. 9 units. Consequently, Rjj j is not a good measure 

of whether any single reaction can be removed and one must try to explain why a quite 

important reaction can be removed without effecting the computed performance. 

One possible explanation is that specific impulse is more directly related to the 

total enthalpy change rate for the whole recombination process than to any single Xtt • 

] 

value. The total enthalpy change rate ^ ^ is the algebraic sum of the individual 

Xtt • values: 

1 


"H, total 


Z 

j=l 


X 




(15) 


where r is the number of reactions in the mechanism. The complete recombination 
process is always exothermic; therefore, Xjj {- ota j will always be negative, since a 
negative Xtt - indicates a net release of energy into the gas stream due to the particu- 
lar reaction. Now it is possible that the mechanism shown in table I uses more reactions 
than are necessary to completely describe the recombination process. Therefore, if an 
important reaction is removed, the X H , values for some of the other reactions may 
become more negative, so that the net effect is to change X H total very little. That 
this is the case for the H + F + M ~ HF + M reaction can be seen from the data in 
table VI. In this table, values of X H j are shown for an assigned chamber pressure of 
100 psia (68. 95 N/sq cm) and O/F of k using three different reaction mechanisms. One 
is the complete set of reactions in table I. Mechanism A is all these reactions except 
reaction I, H + F + M ~ HF + M. Mechahism B is all the reactions of table I except re- 
action Xm, BO + HF ~ BOF + H. The values of X^ tota j are also shown for each 
mechanism. Comparing mechanism A with the complete mechanism shows that the 
X H i values for several reactions (II, III, IV, and X) become significantly more nega- 
tive when reaction I is not present. The total effect gives a value of X ^ which is 

only 1. 5 percent less negative than the value for the complete mechanism. For mecha- 
nism B the removal of BO + HF ~ BOF + H causes Xtt ■ values for reactions XI, XIV, 

H, ] 

XVII, and XXV to become more negative, but causes the values for reactions X and XVIII 
to become more positive. The result of these changes is to make Xjj about 5 per- 

cent less negative than the value for the complete mechanism. Therefore, when the 
magnitude of Xjj to ^- a ^ decreases by 1. 5 percent, specific impulse is negligibly changed 
(less than 0. 5 impulse units). When the magnitude of Xjj decreases by 5 percent, 

specific impulse decreases by about 2 units, a nonnegligible change. This single com- 
parison is not conclusive, but does indicate that the total enthalpy change rate is more 
important than any single Xtt i value in estimating the effect of a mechanism change on 
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TABLE VI. - ENTHALPY EXCHANGE RATE PARAMETERS USING 


DIFFERENT REACTION MECHANISMS 


[Chamber pressure, 100 psia (68.95 N/sq cm); oxidant 
fuel ratio, O/F = 4. ] 


Reaction 

Reaction 

Enthalpy exchange rate parameter, X,, . 

J 

number 







Complete 

Mechanism A 

Mechanism B 



mechanism 



II 

h + h + m = h 2 + m 

-311 300 

-369 800 

-316 800 

I 

H+F + M-HF + M 

-117 800 

0 

-119 100 

IV 

H + OH + M - HgO + M 

-101 900 

-120 700 

-105 600 

xm 

BO + HF ~ BOF + H 

-100 800 

-100 900 

o 

XIV 

BF + OH - BOF + H 

-24 670 

-24 550 

-79 100 

m 

F + H 2 ~ HF + H 

-23 320 

-35 490 

-19 760 

VI 

Hg + OH — H 2 0 + H 

+ 16 200 

+19 140 

+14 250 

XIX 

BF + O - BO + F 

+ 11 370 

+11 240 

-1 726 

xvm 

H + BF 2 - HF + BF 

-8 423 

-8 623 

+3 295 

V 

O+H + M-OH + M 

-7 889 

-9 445 

-8 078 

XI 

BO + F + M ~ BOF + M 

-7 576 

-8 957 

-14 970 

XXV 

BF + 0 2 - BOF + O 

-6 183 

-6 765 

-21 790 

xn 

BF + O + M - BOF + M 

-5 528 

-6 336 

-8 768 

X 

F + OH - HF + O 

+5 283 

-11 130 

+62 000 

XVII 

BO + BF 2 - BOF + BF 

-3 994 

-4 031 

-27 250 

XX 

B0 2 + BF ~ BOF + BO 

-3 399 

-3 528 

-2 384 

vm 

O + H 2 - OH + H 

+2 350 

+ 3 437 

+2 399 

VII 

H + 0 2 - OH + O 

-2 120 

-2 543 

-8 684 

XXI 

BF 2 + HF - BFg + H 

+ 1 946 

+2 029 

+ 1 291 


Sum of X„ .'s 
H, ] 

-687 700 

-677 000 

-650 800 


performance. This idea is consistent with the interpretation of Xjj ^ ^ as 3X1 approx- 
imate indicator of the enthalpy release rate into the gas stream. ^ ^ is then re- 
lated to the kinetic energy of the gas and therefore to specific impulse, which is propor- 
tional to the velocity of the gas stream. 

It should be mentioned that only 19 reactions are listed in both tables IV and VI. 

The additional 6 reactions in table I were used during the calculations, but their contri- 
bution to X H t 0 £ a i an d to computed performance is negligible, so they are omitted from 
the tables. Additional calculations showed that only the first 16 reactions listed in 
tables IV and VI are needed in order to compute kinetic specific impulse within our re- 
quired uncertainty of ±0. 5 pound-second per pound (4. 9 (N- sec) /kg) for all assigned con- 
ditions. Only one of these is really rate controlling. This is the three-body recombina- 
tion H+H + M-^Hg + M. Performance is sensitive to rate constant uncertainties for 


17 


two other important recombinations, H + F + M — HF + M and H + OH + M — HgO + M. 
Most of the remaining reactions are rapid, low energy shuffle reactions which supply the 
free radical species that react in the three- body recombination processes. The most 
important reaction of this type is BO + HF BOF + H, which generates hydrogen atoms. 


Comparison with Hydrocarbon - Fluorine/Oxygen Propellants 

In order to compare the most promising space storable propellant systems, theo- 
retical performance curves are presented for three systems in figure 5. Equilibrium 
and kinetic results obtained in the present work for OFg-diborane are shown along with 
those presented in reference 1 for the methane - Fg/Og and butene- 1 - Fg/Og systems. 
The comparison is made at 100 psi chamber pressure (68. 95 N/sq cm) and a nozzle ex- 
pansion ratio of 40. The fluorine/oxygen mixture giving best performance was used for 
each of the hydrocarbon fuels. It must be noted that the actual nozzle area profiles are 
not the same for the hydrocarbon systems as for the OFg-diborane propellant. This 
should not, however, have a significant effect on the comparison. Figure 5 shows that 
both equilibrium and kinetic curves for the latter system lie well above the correspond- 
ing curves for either hydrocarbon propellant. The kinetically-limited performance level 
for OF 2 -diborane is 13 to 24 pound-seconds per pound (127 to 235 (N-sec)/kg) higher 
than the maximum kinetic specific impulse for methane-fluorine/oxygen. Moreover, the 
actual percentage kinetic loss is less for the OFg-diborane system than for the others. 


4400i — 


4200! 


4000 


3800 


3600 


3400 


Oxygen difluoride-diborane 



Oxidant/fuel ratio, 0/F 


Figure 5. - Theoretical performance of "space storable" propellant systems. 
Chamber pressure, 100 pounds per square inch absolute (68. 95 N/cm^); 
final expansion ratio, 40. (Reference values for all reaction rate con- 
stants. ) 
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The ratio of kinetic to equilibrium specific impulse at the O/F value for peak kinetic 
performance is about 0. 97 for OFg-diborane and 0. 96 for the other two systems. There- 
fore, the OFg-diborane has higher kinetic performance not only because its equilibrium 
performance is higher, but also because its kinetic losses are smaller than for the 
hydrocarbon-fluorine/oxygen systems. This comparison shows the superiority of the 
OFg-diborane system in terms of theoretically computed performance. 


CONCLUDING REMARKS 

The results of the present work have shown that OFg-diborane nonequilibrium per- 
formance is sensitive to only 3 of the 25 possible reactions in the nozzle. Only the rate 
constants for the reactions 


h+h + m-h 2 + m 

H+F + M-HF + M 
H + OH + M ~ H 2 0 + M 

must be known to better than one order of magnitude accuracy to compute performance 
within an uncertainty of ±0. 5 pound-second per pound. Thirteen additional reactions (in- 
cluding those of the boron compounds) must be in the mechanism, but their rate con- 
stants can be in error by one order of magnitude or more without changing specific im- 
pulse more than the stated uncertainty. The hydrogen atom recombination is really the 
one dominant reaction, inasmuch as computed performance is considerably more sensi- 
tive to this rate constant than to any other one. The OFg-diborane system is one of the 
most complicated of the space storables studied so far. It contains more species and 
requires more chemical reactions than do the hydrocarbon-fluorine/oxygen systems. 

The latter systems have at least two dominant reactions (the H + H and H + F recombina- 
tions) whereas OFg-diborane has only one. 

The theoretical results shown here indicate that OFg-diborane has higher specific 
impulse and smaller kinetic losses than the hydro car bon- Fg/Og combinations. Unfortu- 
nately no experimental performance data were found for OFg-diborane which were felt 
adequate for comparison with these computations. Therefore no comparison of experi- 
mental and theoretical results can be made. Uncertainty in the theoretical results can 
be due to other factors beyond any rate constant uncertainties. One of these is the value 
of the heat of formation of oxygen difluoride used to perform the equilibrium computa- 
tions that are the starting point for obtaining the kinetic results. The current value 
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given by the J. A. N. A. F. Thermo chemical tables for the heat of formation of OFg(g) at 
298. 15 K is AHggg = -4. 39 kilocalories per mole (-18.4 kJ/mole). As stated previously 
we have used the value AHggg = 5. 86 kilocalories per mole (24. 52 kJ/mole) recently 
measured experimentally by King and Armstrong (ref. 15). This new value, higher by 
almost 10 kilocalories, increases the theoretical performance level for OFg-diborane ap- 
proximately 7 pound-seconds per pound (69 (N- sec) /kg). Therefore, this combination 
would still be more energetic than the hydrocarbon- Fg/Og combinations even if the value 
AH^gg = -4. 39 kilocalories per mole had been used for OFg. The use of the recently re- 
ported higher heat of formation allows OFg-diborane to compare even more favorably 
with other systems. Of course, this comparison has not considered any of the practical 
problems involved in using either type of propellant combination in a practical engine. 


SUMMARY OF RESULTS 

Theoretical computations have been performed for the oxygen difluoride- diborane 
propellant system to obtain both equilibrium and kinetically-limited performance at 
chamber pressures of 100 and 300 psia (68. 95 and 206. 9 N/sq cm) over an O/F range of 

2. 7 to 4. A study was also made to select the important chemical reactions occurring 
during the nozzle expansion and to evaluate the sensitivity of performance to individual 
variation in the rate constants of these reactions. The results of thie investigation are 
as follows: 

1. Theoretical performance is moderately dependent on the assigned oxidant-fuel 
ratio. For the fuel-rich O/F range used, computed values of vacuum specific impulse 
increase steadily as O/F is increased. The kinetieally limited values for 100 psi cham- 
ber pressure (68. 95 N/sq cm) range from 411 to 422 pound-seconds per pound (4030 to 
4138 (N-sec)/kg). At 300 psia (206.9N/sq cm), rate limited performance varies from 
423 to 434 pound-seconds per pound (4148 to 4256 (N-sec)/kg). 

2. The most important reaction occurring in the nozzle is the hydrogen atom recom- 
bination H+H + M-H^ + M. Increasing or decreasing this reaction's rate constant by 
the estimated uncertainty changes computed specific impulse values by 3 to 6 pound- 
seconds per pound (29 to 58 (N- sec) /kg), or slightly more than 1 percent, over the entire 
O/F range used. 

3. In the low oxidant-fuel ratio range performance is not significantly affected by 
any other reaction rate constant when these constants are changed within the limits of 
their estimated uncertainty. 

4. In the high oxidant-fuel ratio range specific impulse is also affected by uncer- 
tainty in the H + OH + M ~ HgO + M and H + F + M ~ HF + M rate constants. The effect 
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is greater for the first reaction, since the uncertainty is greater for the H + OH recom- 
bination. 

5. The results of the reaction-importance study are in good agreement with the pre- 
dictions of the reaction rating parameter Xtt •, which measures the enthalpy change 

J 

rate in the gas caused by the j reaction. 

6. Theoretical kinetic performance is not sensitive to one order of magnitude un- 
certainty in the rate constants for the reactions involving the boron- containing species, 
although many of these reactions must be present in the mechanism. The most impor- 
tant boron- containing reaction is BO + HF — BOF + H, which generates hydrogen atoms 
for the recombination reactions. 

7. The kinetically-limited performance level for OFg-diborane is 13 to 24 pound- 
seconds per pound (127 to 235 (N-sec)/kg) higher than the maximum vacuum specific im- 
pulse for methane-fluorine/oxygen, the most energetic of the hydrocarbon-fluorine- 
oxygen systems. The superiority of OFg-diborane is due not only to the fact that its 
equilibrium performance is higher, but also to the fact that its kinetic losses are lower 
than those for the methane-fluorine/oxygen system. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 25, 1968, 

129-01-06-01-22. 
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APPENDIX - SYMBOLS 


A 

sd 


( c p1 


c i 

E„ 


h- 


eq 

T k 

Vac 

M 

Mi 

Jl 

P 

R 


R 


r 

T 

V 


* 


W i 


X 


X 


H, total 


AH! 


298 


o 

nozzle crossection area, length 
parameter defined in eq. (8) 
parameter defined in eq. (9) 


heat capacity of mixture based on fixed composition = 
energy/(mass)(K) 

heat capacity of species i, energy/(mass)(K) 


N 


Z 



concentration of species i, mole/volume 
activation energy, cal/mole 
enthalpy of species i, energy/mass 

theoretical equilibrium vacuum specific impulse, (lb-sec) /lb ((N-sec)/kg) 

theoretical rate limited vacuum specific impulse, (lb-sec) /lb ((N-sec)/kg) 

theoretical vacuum specific impulse, (lb-sec) /lb ((N-sec)/kg) 

average molecular weight of mixture 

molecular weight of species i 

Mach number defined in eq. (6) 

pressure, force/area 

universal gas constant, energy/(mole)(K) 


reaction rating ratio = 


X H,i/( X H,j), 


max 


reference length 
temperature, K 
velocity, length/time 

production rate of species i, mass/( volume) (time) 
reaction rating parameter defined by eq. (11) 

net enthalpy change rate for complete nozzle expansion process, defined by 
eq. (15) 

heat of reaction at 298 K for the forward reaction as written in table I, 
kcal/mole (kJ/mole) 
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J.'L 

net conversion rate parameter for j reaction defined by eq. (12) 

net volumetric conversion rate for reaction defined by eq. (12) 

nondimensionalized distance coordinate = distance/r* 
mass fraction of species i 
heat capacity ratio defined by eq. (7) 

stoichiometric coefficients for species i in reaction j; see eq. (13) 
density, mass /volume 
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